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Abstract 

A validity of a recent estimate of an upper limit of charm production in central 
Pb+Pb collisions at 158 A-GeV is critically discussed. Within a simple model we 
study properties of the background subtraction procedure used for an extraction 
of the charm signal from the analysis of dilepton spectra. We demonstrate that 
a production asymmetry between positively and negatively charged background 
muons and a large multiplicity of signal pairs leads to biased results. Therefore the 
applicability of this procedure for the analysis of nucleus-nucleus data should be re- 
considered before final conclusions on the upper limit estimate of charm production 
could be drawn. 
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Measurement of the invariant mass spectra of opposite-sign lepton pairs (dileptons) 
allow to extract information otherwise difficult or even impossible to obtain. Among 
interesting processes which contribute to dilepton production are decays of vector mesons 
(p, w, 4>, J/if), ip), Drell-Yan as well as thermal creation of dileptons, and decays 
of charm hadrons. Decays of pions and kaons are a dominant source of uninteresting 
(background) dileptons which should be subtracted before deconvolution of contributions 
from the interesting (signal) sources is performed. 

Recent analysis of dimuon spectrum measured in central Pb+Pb collisions at 158 
A-GeV by NA50 Collaboration jj]] suggests a significant enhancement of dilepton pro- 
duction in the intermidiate mass region (1.5+2.5 GeV) over the standard sources. The 
primary interpretation attributes this observation to the increased production of open 
charm |lj . In the following theoretical papers other possible sources of the observed effect 
are proposed which do not invoke enhancement of the open charm yield . This suggests 
to interpret the NA50 result as an estimate of the upper limit (about 3 times above pQCD 
predictions) of open charm multiplicity in Pb+Pb collisions at SPS. The above conclusion 
relies, however, on the assumption that the background subtraction procedure used to ex- 
tract signal sources gives unbiased results. In this work we show that this assumption is 
questionable. In particular, an asymmetry in the production of positively and negatively 
charged background dileptons and a high multiplicity of signal pairs lead to the result 
which differs from the one usually assumed in the data interpretation. Our analysis is 
done within a simple model based on the assumptions used to justify the background 
subtraction procedure 0. 

In central Pb+Pb collisions at SPS, due to high multiplicity of produced hadrons, 
the multiplicity of background dileptons is much higher (« 95 %) than the multiplicity 
of signal pairs (~ 5 %). The invariant mass spectra of the Drell-Yan, thermal, and 
open charm contributions are broad and essentially structureless. Consequently their 
extraction requires very precise knowledge of the shape and the absolute normalisation of 
the background distribution. The necessary accuracy can not be reached by calculation of 
the background based on a model. Therefore in order to decrease the systematic error of 
the background estimation a method based on the measured data was developed and used 
in the analysis of dilepton spectra JJ, Q M. In this method the background contribution 

to dilepton spectra is calculated as 2 J (n ++ ) (n ), where (n ++ ) and (n ) are measured 

multiplicities of like-sign lepton pairs. 

The NA50 experiment measured the mean multiplicity of like-sign, (n ++ ) and (n ), 

and opposite-sign, (n + _), muon pairs. One usually distinguishes two classes of muons: 
the "independent" muons coming from decays of pions and kaons (h mesons) and the 
"correlated" muons originating from vector meson decays, Drell-Yan and thermal cre- 
ation of dimuons, and from decays of pairs of charm hadrons. For simplicity of the initial 
considerations let us assume that the correlated muons only come from the decays of 
charm hadrons, which we denote here by D and D. The meaning in which the words 
"independent" and "correlated" used above is the following. Let iV + , N_ be the num- 
bers of positively and negatively charged hadrons (kaons and/or pions) produced in a 
given nucleus-nucleus (A+A) collision. The numbers N + , N_ are independent when the 
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probability to observe them can be factorized: 



P(A + , A_) = P + (N + ) x P_(AL) , (1) 

where P + (N + ) and P_(AL) are the probability distributions for independent observation 
of N + or AL hadrons. Due to charm conservation the numbers of D and D hadrons are 
expected to be equal in each event (AL> = Njy); the production of D and D hadrons is 
correlated. The independence or the correlation of muon sources leads to an independence 
or a correlation of muons originating from these sources. The assumption of approximately 
independent K + and K~ (or tt + and 7r~) production in A+A event is justified by large 
number of different hadron species created in the collision. Then, e.g. the electric charge 
and strangeness of produced K + in a given event could be in fact compensated by many 
different hadron combinations, not just only by K~ . 

Let us denote by cih and old the probabilities that a decay of a single h or D leads to 
a muon inside the NA50 spectrometer. In an event with multiplicities iV + , AL and AL> 
the probabilities to observe n muons of a given sort are binominaly distributed: 

P<«) = „<, -nr, K)" V (i - «*)" + ~ * » (2) 



PM-) = m K) ni (i - a h ) N - - n - . (3) 
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Pc(n c + ) = -j^ c „ (a D r + (1 - a D )»° - ^ , (4) 



Point) = C](N D - , (a D ) n - (1 - a D ) N ° ~ n ~ . (5) 

where n l + , n l _, n c + and nf. are numbers of positively and negatively charged muons from 
"independent" and "correlated" sources. From Eqs. (0-0) one finds 



n\ = ah N + , n % _ = ah N_ , n c + = nl = ajj N D , (6) 
= a h (1 - a h ) N + + a 2 h Nl, (7) 
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a h (1 - a h ) N_ + a\ N 2 



(n%f = (nt) 2 = a D (l- a D ) N D + a 2 D N 2 D . (9) 

We introduce now the probabilities, Ah, Ad, and AhD that muon pairs from, respectively, 
hh, DD and hD decays are detected within the dimuon acceptance. These probabilities 
depend on cuts on the dimuon properties and, for given experimental cuts, on momentum 
spectra of dimuon sources. Assuming that the probabilities A are multiplicity indepen- 
dent, we arrive at the following expressions for the numbers of like-sign and opposite-sign 
muon pairs, for fixed values of iV + , AL and 

^ = A h £ ^ " 1} P«) + A D £ n + {n \ ~ l) P c «) (10) 
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+ AhD E n + n % P i(<) P c«) 

= ~y (i n +) 2 - n +) + ~y (i n %) 2 ~ w+) + A hD n\r}^ 

Ah 2/at2 »r \ , A 
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I (Nl -N + ) + -y- a 2 D (Nl - N D ) + A hD a h a D N + N D , 



n i 



nZT = ^ a£ (iV! - iV_) + ^a 2 D (Nl - N D ) + A hD a h a D N. N D , (11) 



= A E n >- + ^ E n >- p <=«) p <=( n -) ( 12 ) 

+ ^ n\ n c _ PiinX) P c {n c _) + A hD £ n l l < ^(nl) P c (n c + ) 

= A ct\ N + iV_ + A D a 2 D N 2 D + A hD a h a D N D (N+ + AL) . 



Here we have made a simplified assumption that the shape of momentum spectra of h + 
and h~ (as well as D and D) are similar and, therefore, A^ + = A^~ = A^~ = Ah, 
A^d = A^p = A^D = AhD and Aj + = A^ - = A~l~ = Ad (the last equation means 
that possible momentum correlations between D and D are also neglected). Note that 
if there are no cuts on the dimuon properties the above probabilities become equal to 
unity, Ah = AhD — Ad — 1, i.e. assuming all ^-probabilities equal to one in Eqs. (|10[ - 
Ti2"D we count all possible dimuon pairs. However, as soon as one fixes some dimuon 



properties (e.g. an invariant mass of the dimuon pair) all ^-probabilities are evidently 
smaller than unity and their actual numerical values become dependent on the shape of h 
and D momentum spectra and their decay kinematics. Note also that in Eqs. ( |il]| - |T2| ) an 
independence of muon numbers n\ and n l _ is due to assumed in Eq. (|l|) independence of 
N + and iV_ which entered into Pi(n % + ) (|2|) and Pi(n l _) (|3|). A correlation of muon numbers 
n\ and n°_ is due to the correlation of Nd and Njy (Nd = Njj) which entered into P c {n c + ) 
(||) and P c (n c _) (^j) probability distributions. The correlation of n c + and n c _ is of course 
weaker than that for iV^ and Njy, so that n c + are not necessarily equal to n c _ in each event. 

In order to find the final mean multiplicities of the dimuons one should average the 
obtained numbers over all possible values of N + , N_, Nd- To simplify the following cal- 
culations we assume that the relevant multiplicity distributions are Poisson distributions 

N N 

P(N) = — exp(-iV) . (13) 
In this case one gets: 

(n++) £ P(N + ) P(AL) P(N D ) = \ A h a\ i^) (14) 
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(16) 
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Note again that Nd = Njy is assumed in each event and, therefore, there is no independent 
summation over Njy in the above equations. Eqs. (|H}-16) can be rewritten as 
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by introducing the following notations: 
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-a h h 2 _ + - a d D 2 + a m h- D , 
ah h + h- + a d D 2 + a d D + a m D (h A 
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(17) 

(18) 
(19) 
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Parameters a d and a m are therefore the probabilities to observe two muons from the 
corresponding hadron sources (these probabilities are a^, a 2 D and othOtD) within experi- 
mental cuts on muon pair properties (these cuts lead to additional factors Ah, Ad and 
A hD ). 

In the experimental procedure the background contribution to the dimuon spectrum 
is calculated as: 



{nl 9 -) = 2 yj(n ++ ) (n__> . 
The number of signal (/z + , /i~)-pairs is assumed to be: 
(n s + 9 l) = (n+_) - (njf ) . 



(22) 



(23) 



It is expected that the subtraction procedure pp) cancels out all false (/i + , /i~)-pairs i.e. 
the pairs from hh and hD decays, and that (n + 9 _) is proportional to the multiplicity of D 
hadrons: 



(nl 9 l) = a d D . 



(24) 



Let us consider some properties of the subtraction procedure (|23| ) by discussing two 
simple examples within the model. 
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Example 1: We assume that there is no contribution from D-decays. In our model 
this assumption can be introduced by setting «d = 0. Consequently a d = a m = and 
Eqs. (|TT| - |T9|) result in: 



1 1 

(n ++ ) = - a h h\ , (n+_) = - a h h 2 _ , (n + _) = a h h + h_ . (25) 

Using Eq. ( [23D one obtains that (n + 9 l) = 0, i.e. the measured signal multiplicity is equal 
to zero as expected in the case of absence of dimuons from the correlated source. This 
result is valid for any value of h + and 

Example 2: In this example we assume that there are correlated dimuons a d D > but 
the number of positively and negatively charged background hadrons is equal (h + = h_ = 
h). Under these conditions Eqs. ( |l7) - |T9l) can be rewritten as 

(n ++ ) = = - a h h 2 + - a d D 2 + a m h D , (26) 

(n+_) = a h h 2 + a d D 2 + a d D + 2 a m h D . (27) 

Eq. ( p3|) gives (ti+L) = a d D which agrees exactly with the expectation (f24|). 

Finally we consider the general case, i.e. a d D > and h + ^ h-. This last condi- 
tion corresponds to the relation between pion and kaon average multiplicities measured 
in heavy ion collisions: (tt~) > (tt + ) and (K + ) > (K~). From Eqs. ( |T7| - |2"3D by straight- 
forward calculations one finds 

(n S + 9 l) = (n+_) - y/((n + _) - a d D) 2 + 7 D* , (28) 

where 

7 = (a h a d - a 2 n )j (h+ - h_f 



It is easy to see that for an = and/or h + = h_ = h one gets 7 = 0, and the results 
obtained in Examples 1 and 2 are reproduced. We repeat that in the absence of cuts on 
the dimuon properties one has Ah = A^d = Ad = 1. Therefore, a>h a d — = (i.e. 
7 = 0) and consequently we have again unbiased estimate of the mean multiplicity of D 
mesons. In general, however, the result differs from the expected one (f24j). The presence 
of experimental cuts on dimuons (e.g. one fixes the dimuon invariant mass in the region 
M = 1.5 -T- 2.5 GeV) causes that the probabilities Ah, AhD and An are smaller than unity, 
destroys the equality A^Ad = A\ D and, therefore, leads to non-zero value of 7. With cuts 
on dimuon properties the experimental number of signal pairs {n+L) is not equal to a d D. 
By fitting a d D* to (n s + 9 l) one finds the spurious number of D hadrons which we denoted 
by D*. There are two distinct cases. 

Case 1: an a d — o? m < , (7 < 0) . 

The experimentally measured signal, (n+1) fl2"Bp, is larger than the expected value a d D 
and therefore the extracted spurious number of D hadrons is larger than the true one 
(D < D*). 

Case 2: a h a d — a 2 m > , (7 > 0) . 

The experimentally measured signal, (n+ 9 l) (|28|), is smaller than the expected value a d D 
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and therefore the extracted spurious number of D hadrons is smaller than the true one 
(D > D*). 

In the NA50 analysis of the dimuon spectra in terms of the open charm enhancement 
the used background subtraction procedure was checked for two different cases. First of 
all, it was shown to work correctly for simulated central Pb+Pb collisions at 158 A-GeV. 
However in this simulation correlated (signal) muon sources were not included. Thus this 
check is equivalent to our Example 1, for which the procedure works exactly. Secondly 
the open charm yield was extracted for p+A interactions and it was shown to agree with 
the yield from direct measurements. Eq. (^) and Example 1 show that the deviation 
from the expected result decreases with decreasing multiplicity of D hadrons. Thus the 
success of the procedure applied to p+A interactions does not proof its applicability to 
Pb+Pb collisions in which multiplicity of D hadrons may be higher even by a factor of 
about 10 4 H, |. 

Note that our results are obtained in a highly simplified model. The assumptions 
concerning independent production of background muons (Eq. ([[])), the Poissonian mul- 
tiplicity distributions of hadrons (Eq. (|T3|)) and the absence of D meson momentum 
correlations {Ap + = = = Ap) seem to be questionable or even incorrect. Dis- 
cussion of the possible additional biases introduced by these effects is beyond the scope 
of this paper. We also do not attempt here to calculate numerical values of A h , A hD , A D 
for the specific NA50 experimental acceptance of dimuon pairs. 

We close the paper by concluding that the applicability of the background subtraction 
procedure widely used in the analysis of dilepton spectra in nucleus-nucleus collisions 
should be reconsidered. In particular final statement on the upper limit of the open 
charm multiplicity in central Pb+Pb collisions at 158 A-GeV resulting from the analysis 
of the dimuon spectrum requiers further studies in order to quantify a magnitude of the 
bias. They should include numerical simulations of the specific experimental set-up and 
consider various particle production models. 
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